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Abstract. Some physical imprints of quintessence scalar fields on dark matter (DM) clustering are illustrated, and a compar- 
ison with the concordance model ACDM is highlighted. First, we estimate the cosmological parameters for two quintessence 
models, based on scalar fields rolling down the Ratra-Peebles or Sugra potential, by a statistical analysis of the Hubble dia- 
gram of type la supernovae. Then, the effect of these realistic dark energy models on large-scale DM clustering is established 
through N-body simulations. Various effects like large-scale distribution of DM, cluster mass function and halos internal ve- 
locities are illustrated. It is found that realistic dark energy models lead to quite different DM clustering, due to a combination 
of the variation of the equation of state and differences in the cosmological parameters, even at z = 0. This conclusion contra- 
dicts other works in the recent litterature and the importance of considering more realistic models in studying the impact of 
quintessence on structure formation is highlighted. 



INTRODUCTION 

In recent years, there have been increasing evidences in favour of an unexpected energy component — often called 
dark energy — which dominates the present universe and affects the recent cosmic expansion. The first evidence 
came from the Hubble diagrams of type la supernovae [1, 2, 3, 4] and were thereafter confirmed by the measurements 
of Cosmic Microwave Background (CMB) anisotropics (cf. [5] and the other companion papers), the large-scale 
distribution of galaxies [6, 7] and indirectly by the observed properties of galaxy clusters [8] and peculiar velocities 
[9]. This acceleration can be explained by the present domination of an exotic type of matter — often called dark 
energy and being characterised by a negative pressure — over all other ordinary types of energy with positive equation 
of state (eos) (0 = p/p. 

However, several questions have arised on the nature of the underlying physical process behind this accelera- 
tion. The usual and most simple explanation is that the acceleration is due to a positive cosmological constant, which 
acts like a fluid with a constant eos Oa = ^1- This interpretation leads nevertheless to a few intricate problems. 
First, the value of the cosmological constant has to be extremely fine-tuned to account for the observed amount of 
accelerating energy: A w 5 x IQ-^'^GeV^ for Q.Afi = 0.7 (in natural units: /i = c = 1). We refer the reader to [10] 
for a review on this fine-tuning problem and to [11] for an interesting alternate interpretation. The situation about 
the cosmological constant has been made even worst by the discovery of dark energy in recent years because the 
observational evidences claims for a value of vacuum energy density pA of the order of the present critical density 
of the universe. In addition to the fine-tuning problem evoked above, the interpretation of dark energy in terms of 
a cosmological constant coming from non-vanishing vacuum energy also faces this coincidence problem to explain 
this observed value of A. These problems concerning the cosmological constant have made the need of different 
interpretation of dark energy more and more crucial. 

A negative eos co — p/p can be easily obtained in the framework of scalar fields coupled to gravity. This sce- 
nario, often called quintessence, relies on the non-conformally invariant dynamics of such a neutral scalar field 
coupled to gravity. An impressive number of such quintessence scenarii have been suggested in the past few years, by 
Ratra & Peebles [12], Brax & Martin [13], Barreiro, Copeland & Nunes [14], Frieman et al. 1995 [15] amongst many 
others. An accelerating phase can be easily achieved if the potential energy of the field dominates its kinetic energy, 
leading to a negative pressure. Even better, an appropriate choice of the self-interaction potential for the scalar field 



yields to the existence of an attractor which fortunately rules out the fine-tuning problem. This is known as "tracking 
potentials" (see [16]). But the price to be paid for it is the reliance on a hypothetic scalar field ruled by a particular 
self-interaction potential. However, although some physical motivations have been presented for such models, based 
on supergravity, compactified extra- dimensions or spontaneous symmetry breaking and so on, a definitive, physical, 
interpretation has still to be provided. 

Due to those various possibilities, it appears crucial to find a way to trace back one particular model, for exam- 
ple a given potential, from observations. This question is quite intricate as it has been shown that present (and even 
future) data coming from high-redshift supemovae observations allow too much degeneracy between very different 
eos (see [17, 18]). Data coming from supernovae cannot therefore discriminate alone different quintessence scenarii 
and even different analytical approximations of the evolution of the eos fi)(z), at such low redshifts (z < 1.5). Such an 
analysis has therefore to be completed by including physical informations at higher redshifts in order to enlighten the 
dark energy's obscure details. In [19], the authors have performed an exhaustive study of CMB anisotropics within 
the framework of two well-known tracking quintessential models with the Ratra-Peebles [12] (RP) and Sugra [13] 
potentials. They found that such models are fully compatible with BOOMERanG and MAXIMA- 1 data. Indeed, the 
dark energy was almost negligible at the recombination epoch (z ~ 10^) so that most of its effects are concentrated 
in post-recombination anisotropies. This yields to slight effects in the CMB anisotropies, which are finally almost 
insensitive to a dependance on the redshift of the eos (0{z) at recent epochs. 

Large-scale structures formation has been considered with a growing interest in recent years for the dark en- 
ergy debate because of intermediate redshifts (0 < z < 1000) this process covers. As early as in 1999, Ma and his 
collaborators [20] performed N-body simulations to determine the mass power spectrum with an effective dark energy 
with a constant eos throughout the cosmic expansion. Other works by Bode [21], Lokas [22], Munshi et al. [23] 
concentrated on cluster mass functions but all proceeded with the same, rather restrictive, assumption of a constant 
eos. However, analytical works on the linear and second order growth rate of the fluctuations as well as the shape 
of non-linear power spectrum using analytical approximations were performed by Benabed & Bernardeau in [24] 
with RP and Sugra quintessence models. More recently, Dolag and his collaborators [25] studied numerically the 
concentration parameters of dark matter halos and their related evolution in the framework of several dark energy 
models with both constant and evolving eos (RP and Sugra models). Klypin and his collaborators [26] also performed 
N-body simulations on both constant and evolving eos (RP and Sugra). The last found slight differences in the mass 
power spectrum and the cluster mass function at z = that render them almost undistinguishable from the predictions 
of a standard ACDM model. They also discovered more significative differences at higher redshifts in those quantities 
as well as different halo profiles at all redshifts. However, they based their analysis solely on models with same 
cosmological parameters and in poor agreement with supemovae data as we will show thereafter. In [27], the authors 
used the same cosmological models as in [26] to study galaxy distribution as a test of varying eos. They found again 
large differences at large redshift between the dark energy models which are reduced substantially by geometrical 
factors. They concluded that the galaxy redshift distribution should be used as a test of varying eos instead of the 
cluster redshift distribution which leads to too slight differences. However, we wiU show here that considering realistic 
quintessence models — with different cosmological parameters selected by a Hubble diagram analysis — enlarges 
the discrepancies due to the different cosmic evolutions, even at z = 0. 



According to the cosmological principle, the large-scale Universe is roughly homogeneous and isotropic and can thus 
be well described by the Friedmann-Lemaitre-Robertson- Walker (FLRW) line element: 



where a{t) is the scale factor, k is the signature of the curvature and dO. is the solid angle element. We will make use 
throughout this paper of the Planck units system in which h = c=l, mpi = 1/ VG = 1.221 1 x lO'^GeV . The Universe 
is assumed to be fiUed with different matter fluids — pressureless matter (composed of ordinary baryonic matter and 
cold dark matter), relativistic matter (photons and neutrinos) and the mysterious fifth component, caUed quintessence. 



ACCELERATED COSIMOLOGICAL IMODELS 




(1) 



that is hoped to account for the observed cosmic acceleration. The expansion rate of the Universe, characterized by the 
Hubble parameter H =d/a,\s dictated by the Friedmann equation: 
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where the subscript zero stands for nowadays value and where the parameters = pi/Pc represent the density of the 
different components of the matter fluids, expressed in critical units. More precisely, the subscripts m, r, K, A and 
Q indicate respectively the contributions of the pressureless matter, radiation, curvature, cosmological constant and 
quintessence to the cosmological expansion. In equation (2), pg stands for the density of the quintessence fluid for 
which there is, in general, no analytical expression in terms of the scale factor a. For the sake of completeness, we also 
remind the reader about the acceleration equation: 

2^ = -Hi (£l^^,^+2n,^,<) (1+30)2), (3) 

where (x)q is the ratio between the energy density pg and pressure pq of the quintessence fluid. 

We will assume here that the quintessence fluid is actually constituted by a neutral (real) scalar field Q{t) which 
couples only to ordinary matter through its gravitational influence (minimally coupled scalar field). The dynamics of 
such field is ruled by the Klein-Gordon equation written in the metric (1): 
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where <j) is the field expressed in units of the Planck mass iQ{t) = mpi(j){t)), and V{(j)) is the field's self-interaction 
potential. The couphng of the quintessence field to gravity is achieved through its energy density and pressure: 

2 

PQ = ^f + V{,j)), 

2 

fix 

PQ = (OqPq = ^<I>^-V {(!>)■ 



At this stage, the coupled dynamics of gravitation and the quintessence field is only determined by the choice of 
a particular potential and initial conditions, i.e. a couple (g,, Qi), in the early universe, for example at the end of 
inflation (a ~ 10^^^). Convenient choices of self-interactions are those of so-called tracking potentials (see [16] 
and references therein) that account for the observed cosmological parameters today from a huge range of initial 
conditions, spread over dozens of order of magnitudes, in the field phase space. 

Two well-known examples of such potentials are the Ratra-Peebles (RP) inverse power law [12], 

= (5) 

where a > and X are free parameters, and the same potential affected by a radiative correction in supergravity due 
to Brax& Martin [13]: 

ys{^) = ^e''^\ (6) 

which we will refer further to the Sugra model. The RP model was originally suggested to mimic a time-varying 
cosmological constant, which is actually accomplished when the scalar field is frozen by the expansion at a high 
energy scale, typically the Planck mass, where the potential is nearly flat. However, as Brax & Martin noticed in [13], 
this state of the field should correspond to a supergravity regime {Q{t) « mpi) and the potential has therefore to be 
corrected accordingly. 



The particular shape in inverse power law of these potentials is hoped to be justified by energy physics (see for 
example [13] and references therein). In this scenario of quintessence, the field energy density remains subdominant 
in most of the cosmic evolution, especially during the radiation dominated era, until we reach the coincidence epoch, 
which depends exclusively on the potential shape, due to the tracking property. 

Let us now estimate the cosmological parameters Q.m,o and Q.Qfi to be used for quintessence models by analysing 
Hubble diagrams of a type la supemovae data set from the SNLS collaboration [4]. The moduli distance versus 
redshift relation of standard candles is given by 

H{z)=m-M = 5logiQdL{z), (7) 
where diiz) is the luminous distance (in Mpc) given by 

ddz) = (1 + z)Ho^'^l sin (^4,0 £ ^) ' ^^.o > « 
dz' 

= {i+z)Ho —-■,if^K.o=0 (8) 
Jo E{z') 

= {l+z)Ho\£lKfl\~^ sinh (^\^K,o\ ^ 1^ ^) Uf ^Kfl < 

(9) 

for a closed, flat and open FLRW background respectively [28] iE{z) = H{z)/Hq is the dimensionless Hubble 

parameter). In order to avoid considering Hq as an additional parameter in the statistical analysis, we marginalize 
the estimator with respect to this parameter (see also [17]). We will therefore use the following estimator [17] 

^^=A-^+log(|) (10) 

where A = (^'"hAzi)-^.^'th{^))\ q = ^/v^^ ^jjcJ C = ijLi ^ with N the number of data, cr,- is the 

uncertainty related to the measurement Hobsizi) and is the theoretical prediction for the distance moduli. We 
choose to fix the value of the parameter a so that Q,mfl and Q,Qfl will be the only free parameters^ The energy scale A 
in the quintessence potential is here determined so that the desired couple (0,^,0,^0,0) is retrieved. 

Table 1 gives the values of X'^ for the best fit RP, Sugra and hCDM models as well as the evaluation of this 
estimator for the quintessence models used in [26, 27]. These last models are characterized by an energy scale of 
IQ^GeV of the quintessence potential (which corresponds to a fa 4) and the same cosmological parameters. Their 
agreement with SNLS data is very poor as they he well outside the 2 — (7 confidence region. Figure 1 represents 
the cosmological evolution of the models in Table 1. The age of the Universe lies between 13.9 (ACDM model) 
and 15.9 (RP) billion years. Also represented are the models previously used in [26, 27] to study the impact of 
quintessence on structure formation (see also Table 1), which are quite different of the realistic dark energy models 
used here. Figure 1 also illustrates the cosmological evolution of the eos (Oq = Pq/Pq for the best fit RP and Sugra 
models. This shows that the same adequacy to data of a Hubble diagram at low redshifts z < 1 — 2 is obtained with 
very different the eos, leaving this quantity very difficult to determine from supemovae data alone as already shown 
in [17]. There is a degeneracy in the analysis of Hubble diagram as the same cosmic acceleration can be obtained 
either with a small quantity dark energy (£1q small) providing strong acceleration ((Oq <^ 0) (for example, this is 
the case of the cosmological constant here) or a larger amount of dark energy exerting weaker acceleration (this is 
the case of quintessence models here for which (Oq > —1). Table 1 gives the coincidence and acceleration redshifts 
Zc iPgizc) = Pm{zc)) and Za respectively for the best fit models of Table 1, iUustrating that quintessence models 
dominate and accelerate earlier than a cosmological constant. As their accelerating power is smaUer {(Oq > — 1) 



In the case of the Sugra model, this is well-motivated by the fact that the exponential correction in equation (6) is dominant in the potential at 
late redshifts (when (p w mpi), leaving the Hubble diagram analysis insensitive to a variation in the parameter a. This was proved a posteriori by 
Caresia and his collaborators (Caresia, Matarrese & Moscardini, 2004). 



TABLE 1. Cosmological parameters and values of on the 
considered data set for the different models. The coincidence and 
acceleration redshifts Zc andza as well as the starting redshift Zin 
of the N-body simulations are also given 
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RP a = 0.5 


0.23 


0.91 
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0.71 


0.99 


33 


Sugra a = 6 


0.2 


0.84 


116 
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41 


A (SNLS) 


0.27 


0.75 
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0.41 


0.78 


37 
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0.3 


0.7 


168 


0.82 


0.08 




Sugra* 


0.3 


0.7 


127 


0.43 


0.51 
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FIGURE 1. Cosmological evolutions of the scale factor (left), acceleration factor (center) q = da/d^ and eos (Oq = Pq/Pq for 
the models in table 1. The toy models used in [26, 27] are given in black curves 



than a cosmological constant, they need more time to provide the same cosmic acceleration observed in supernovae 
data. Therefore, the statistical analysis of Hubble diagram in this section leads to a prediction of closed cosmologies 
(Q.tot = + > 1) for quintessence (see table 1). This conclusion makes the cross-analysis with CMB data very 
interesting for quintessence as this could rule out models with too large Cltot, but we will leave this point for further 
studies. Figure 2 represents the 68 and 95% confidence level contours for the RP and Sugra quintessence models (a 
similar figure for the concordance model ACDM can be found in [4]). The statistical analysis has been performed 
with only two degrees of freedom, the density parameters £2„, and £2g, as the power-law a in (5) and (6) has been 
fixed. In the Sugra model, we took a = 6 so that the energy scale X of the potential is of the order of lO^GeV, a scale 
related to high-energy physics. In this model, the exponential correction in (6) flattens the potential which allows a 
good adequacy to data. With such an energy scale A of lO^GeV (a « 6), the RP potential would not be flat enough to 
provide the required cosmic acceleration and we reduced a to 0.5 for the RP model to make the model compatible 
with data (A « leV). We find Q.,„ < 0.38 (il,,, < 0.31) and Q.q < 1.44 (Hg < 1.21) at the 95% confidence level for the 
RP (Sugra) model. 

To conclude this section, we state that the statistical analysis of the Hubble diagram of type la supernovae favours 
closed quintessential universes. In general, more dark energy is needed to provide the same cosmic acceleration with 
quintessence scalar field than with a cosmological constant. Furthermore, assuming the same cosmological parameters 
for quintessence and concordance model will lead to a bad approximation of the recent cosmic expansion (see table 
1 and Figure 1). That is why a preliminary analysis is unavoidable to determine the best density parameters before 
performing numerical simulations of structure formation that are expected to provide additional constraints on dark 
energy. We will now use the first three models of table 1 to evaluate the discrepancies between cosmological constant 
and quintessence on structure formation in dark energy-driven universes. 
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FIGURE 2. 68 (red) and 95% (black) confidence contours for the RP (left) and Sugra (right) models. 



IMPRINTS OF QUINTESSENCE ON LARGE-SCALE DARK MATTER CLUSTERING 

Our aim is to give an overview of some impacts of quintessence on large-scale structure formation, mainly on cluster 
properties like mass functions and internal velocity dispersion. We will treat here realistic quintessence models that 
account for the Hubble diagram of type la supernovae while other works like [26, 25] made use of toy models 
for quintessence^ and focused mainly on the cluster mass functions. We will study dark matter collapse on large- 
scales through N-body simulations (Particle-Mesh algorithm). We will consider 256^ particles moving into a box of 
comoving length with 256^ grid cells (h = 0.7), with the baryon density parameter equal to 4%. At 
scales lower than IQQh^^Mpc, it has been shown in [19] that a ACDM power spectrum can be used for the initial 
conditions (the shape of the spectrum in a ACDM and in RP and Sugra quintessence models is very similar at those 
scales). We use the WMAP3 normalisation [29] for each of the simulation, with nns density fluctuation level within 
spheres of Sh^^Mpc radius = 0.74. The initial redshifts (given in table 1) of the simulations are determined such 
that the linear evolution of the filtered dispersion dg reaches today this value. The hnear evolution of the matter density 
contrast is given by the solution of the following equation (see [28]): 

d^5"+(d + 2-]8'--^ = Q, (11) 
\ a J 

where we assumed oq = I, A-Ttmp^Pm o — 1 and where a prime and a dot denote a derivative w.rt. the scale factor a 
and time f, respectively. The analytical solution under integral form, originally derived by Heath in 1977 for universes 
with matter and cosmological constant, is no longer valid in the case of quintessence, even with a simple constant 
eos. We have therefore to solve numerically equation (11) with initial conditions given by requiring that the density 
contrast evolves linearly with the scale factor a after recombination z ~ 10^ (matter-dominated universe). 

The left plot of Figure 3 represents a slice of the dark matter 3D density field poM for the best fit ACDM 

model at z = 0. The other plots in Figure 3 (center and right) illustrate the absolute differences Pum"^ / PmF'^ ~ ^ ^'^^ 
the same slice with respect to the ACDM model for the Sugra and RP models respectively. The bright (dark) regions 
in those plots indicate where the quintessence models are more (less) structured than the ACDM. We can see that the 
ACDM is more clustered in the filaments than the quintessence models (dark regions in the central and right parts of 
Figure 3) leaving the last more structured outside the filaments (bright regions). Let us now examine the discrepancies 
in the effect of quintessence or cosmological constant on dark matter halos. Here, they are located thanks to the 
Friend-Of-Friend algorithm, which groups inside a same cluster all the particles that are within a certain fraction of 
the mean distance between particles from each other We will consider clusters of at least 100 particles and a reference 
distance used to identify the cluster of 20% of the mean distance between particles. Figure 4 (left) represents the 
number density of clusters with a mass above a given value Mq at z = and z = 1 for the ACDM models, while the 
central and right part of Figure 4 illustrate the ratio between the number density in models with quintessence and 
the number density of ACDM. The differences are of the order of 50% for high masses, even at z = 0. In the toy 



They do not account for supernovae data. 



FIGURE 3. Dark matter density field in a slice of the ACDM model (left) and absolute differences for the Sugra (center) and RP 
(right) models. The color bar indicate the degree of discrepancy with a ACDM model in percent. 




FIGURE 4. Mass function of dark matter halos in the ACDM model at z = (straight line) and z = 1 (dash-dotted line) (left). 
The mass functions for the quintessence models at z = (center) and z = 1 (right) are illustrated in units of the mass function of the 
ACDM model. 



models used in [26], the mass functions at z = were found almost undistinguishable. Therefore, using the right 
cosmological parameters for quintessence models, as suggested by distance-redshifts measurements, increases the 
differences on the mass functions at z = found in [26] for toy models. Another important difference with [26] is the 
following: the authors also found that the RP model was more clustered than the Sugra which was more clustered 
than the ACDM. This order is, however, due to their assumption of same cosmological parameters which yields that 
the differences in structure formation are only due to different cosmic expansion. This order reflects the acceleration 
provided by quintessence : their RP model accelerates less {(Orp ~ —0.4) than their Sugra model ((Osugra ~ —0.8) 
and the cosmological constant ACDM (coa = —1). Here, we find that the ACDM model is more clustered than the 
RP and Sugra quintessence models, and this fact is due to different cosmic expansion and the different cosmological 
parameters required by distance-redshifts measurements. 

Let us now examine the internal velocity dispersion of dark matter halos, given by a?, = 1 /NYJiLi 1 1"*^; ~ "i^cmI P 
(with the number of particles in a given cluster and ~V^cm the center-of-mass velocity). This is done in Figure 5: 
on the left the number density of of halos with velocity dispersion higher than a given value is given for the ACDM 
model at z = and z = 1 . On the center and right part of the figure, the number densities for quintessence models 
are given in units of the ACDM amount. Quintessence lead to higher internal velocity dispersion than the ACDM at 
z = 0. This can be justified by the fact that these models produce more lighter halos at z = which are less virialised. 
At z = 1, the differences with a ACDM are smaller, which means that more light halos are formed in quintessence 
models between z = and z = 1 (see also Figure 4). 
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FIGURE 5. Number density of halos with internal velocity dispersion larger than a given value vq in the ACDM model at z = 
(straight line) and z = 1 (dash-dotted line) (left). The same number densities for the quintessence models at z = (center) and z = 1 
(right) are illustrated in units of the number density of the ACDM model. 



CONCLUSIONS 

In this paper, we have illustrated some impacts of quintessence on the distribution of dark matter on large-scales. 
We have shown through to an analysis of Hubble diagram of type la supernovae that the cosmological parameters 
(ilm OjfiQ.o) for quintessence models are quite different than in the concordance ACDM model, as might be expected 
from the different cosmic acceleration quintessence models provide. Cosmological models with quintessence are 
spatially closed and exhibit a slightly smaller amount of matter than the concordance model with cosmological 
constant. The effect of the necessary different cosmological parameters was underestimated in the litterature on the 
imprints of quintessence in structure formation. We have shown that this lead to more important deviations between 
the predictions of a ACDM and quintessence models at z = 0. As well, considering the same cosmological parameters 
for quintessence and cosmological constant does not fit Hubble diagram data but it also introduces a bias in the 
interpretation of the imprints of quintessence on the large-scale distribution of dark matter With the same cosmological 
parameters, quintessence models were found in [26] to produce more structures and to be more clustered than a ACDM 
scenario. However, once the correct cosmological parameters are taken into account, we have shown that this is exactly 
the opposite. In this paper, we have shown that quintessence yields to more structures outside the filaments, lighter 
halos with higher internal velocity dispersion. The present study illustrates how dark matter is sensitive to an expansion 
driven by quintessence or cosmological constant. These differences on the dark matter collapse, even at z = 0, will 
lead to observable imprints on galaxy formation. An interesting point for further works will be to establish constraints 
on the nature of dark energy from galaxy properties through numerical simulations. 
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